Dynamical model of steadily forced isotropic turbulence 
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Abstract 

A dynamical model is proposed for isotropic turbulence driven by steady forcing that yields a vis- 
cosity independent dynamics for the small-scale (inertial) regime. This reproduces the Kolmogorov 
spectrum for the two-point velocity correlation function in the fully developed (stationary) stage, 
while predicting intermittency corrections for higher order moments. The model also yields a tran- 



sient stage with a power-law time evolution. The crossover time to fully developed turbulence 

scales with the turbulent system size as ~ L ' s . The physical origin of the transient behavior is 
explained. 
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I. INTRODUCTION 



With regard to the problem of turbulence, the nonlinear inertial term in the Navier-Stokes 
equation is responsible for energy redistribution among the constituent length scales, as well 
as for the so-called sweeping effect, which couples any given scale with all the larger scales 
without redistribution of energy. Attempts to understand these nonlinear effects have led to 
different approaches. For instance, the nonperturbative nature of the sweeping interactions, 
which has been shown to mask the classical scale-invariant theory of Kolmogorov Q], 
has led to the exploration of nonperturbative alternatives (e.g. 3, Q 0, 0, El)- In forced 
systems, other alternatives such as the Burgers equation (e.g. [8|, |3, [l0|), shell models (e.g. 



11 1) or their variants |l2| have been also explored. Some attention has been paid, as well, to 



the possible depletion of nonlinearity in parts of the real space ^J. Yet another approach, 
based on the idea of inventing an effective model for the nonlinearity, has in turn led to 
a large collection of phenomenological models Q]. The model we shall introduce here for 
isotropic turbulence falls more or less into the last category. 

Because isotropy everywhere implies homogeneity, by isotropic turbulence we mean, of 
course, homogeneous isotropic turbulence. In isotropic turbulence, in contrast with shear- 
flow turbulence, no mean shear stress can occur and consequently the mean motion of the 
fluid is uniform. Isotropy is a statistical property of the fluctuating part, u(x, t), of the 
fluid's velocity field, i.e., a property of the velocity field with respect to an observer moving 
with the mean motion of the fluid. The isotropic field u(x, t) satisfies 

< u(x, t) >= 0, < Mj(x, t)Uj(x., t) >= - < U 2 (x, t) > 5ij. 
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Assuming the flow incompressible, in the reference frame of the co-moving observer, the 
Navier-Stokes equation takes the form 

d t u+ (u.V)u = z/V 2 u- -Vp + f, (1) 

where p(x, t) is the fluctuating part of the pressure field (< p >= 0) and f(x, t) represents 
the fluctuating stirring force (< f >= 0). The pressure term is in effect a nonlinear term; 
it can be eliminated via the incompressibility condition in the standard manner (e.g. j3j), 
at the price of introducing additional (and more complicated) nonlinearity in the equation. 
In view of the isotropy, it suffices to consider the dynamics of the velocity profile in a fixed 
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direction. Thus, let v(x, t) and r?(x, t), respectively, denote the fluctuating velocity and force 
components in an arbitrary direction. Equation JIJ in essence reads: 

d t v = — V.jvisc + n + nonlinear terms, 

where j V i SC = —vVv is the viscous current and represents the momentum flux due to viscous 
stresses. The viscous current acts to oppose the spatial variations of v and, thus, has a 
smoothing effect on its fluctuation-induced irregularities. Our approach is based on the idea 
of inventing an effective model for the nonlinearity by augmenting j„ isc by a suitable 'current' 
term j fw .f, according to 

8 t V = -V.(.W C + }turb) + V- (2) 

Of course, j tU rb must satisfy < j t urb >= (up to a trivial constant). Physically, the augment- 
ing inertial current jturb is supposed to represent the momentum flux due to the collective 
action of turbulent fluctuations. The competition between the two currents j V i SC and i tU rb 
will be central to the turbulence problem. In equation © the fluctuating fields n and } tU rb 
are to be deduced from physical arguments and assumptions, rendering the basic equation 
of the model. The model we so present is new and will be shown to yield, through the 
competition of the two currents, a small-scale regime in which the dynamics is independent 
of viscosity. This regime corresponds to the inertial regime of the steadily forced isotropic 
turbulence. The equation associated with the inertial regime can be treated analytically by 
the dynamical renormalization group method to yield the scaling behavior of the two-point 
velocity correlation function. This reproduces the celebrated Kolmogorov scaling law [2| for 
the fully developed (stationary) stage, while yielding the power-law time evolution ~ t 2 ' 11 
for the early transient stage. Furthermore, the model predicts intermittency behavior for 
higher order correlation functions in the fully developed stage. The crossover time from the 
transient to the stationary behavior scales as ~ L 11 / 3 , where L is the turbulent system size. 
The physical origin of the transient behavior will be explained. 

Interestingly, our model equation for the inertial regime (equation (|SJ)) also appears in 
the context of surface growth phenomena j3] as a certain growth equation (see section III). 
While the similarity between the characteristics of surface growth models jl^i Il7|. in par- 
ticular of this growth model Q], and turbulence has been previously noticed, this work 
presents a first instance of derivation of such a model from pure dynamical considerations. 
This provides another perspective from which the model can be viewed. The fact that the 



Kolmogorov spectrum for isotropic turbulence has been derived more or less directly from 
the Navier-Stokes equation (e.g. P, 2l3|), does not hinder the importance of other relevant 
model equations; after all the whole machinery of the Navier-Stokes equation may not be 
required as far as the inertial range dynamics is concerned. In fact, as is well known, in the 
zero viscosity (infinite Reynolds number) limit which is relevant to inertial range dynamics, 
the Navier-Stokes equation shows scaling symmetry for the velocity profile. This motivates 
us to put forward an analogy between interface growth phenomena and the inertial range 
dynamics. The analogy is based on the observation that the turbulent velocity profile may be 
viewed as an evolving interface that attains a self-affine form due to the roughness (width) 
generated by the collective action of turbulent fluctuations. The momentum flux due to 
these fluctuations, as shown in section II, can be obtained from dynamical considerations 
yielding the model equation. 



II. THE MODEL 

In view of the degree of control over the stirring forces exerted in an actual turbulent 
flow, it seems plausible to assume that in microscopic scale the forcing is not correlated in 
space and time, rj may thus be represented by an uncorrelated noise term, with the standard 
properties 

< r/(x, t) >= 0, < r/(x, t)r](x.', t') >= 2D 5 3 (x - x.')5{t - t'). 
Equation ©, then, reads: 

d t v = -V.j + 77, 

where j = ] V isc+'hurb is the total current. At every point v and its derivatives are independent 
variables that can develop (independent) irregularities through fluctuations and the total 
current acts as to oppose the resulting (independent) spatial variations by redistributing 
the irregularities. The viscous component j« sc of the current, as mentioned before, serves 
to iron out the spatial variations of the velocity field v itself. The rest (irregularities of the 
derivative fields, most notably the energy dissipation rate) is to be compensated for by the 
effect of the inertial component jturb- We, therefore, hypothesize that the turbulent current 
has a similar smoothing effect on the spatial variations of the derivative fields, most notably 
the velocity gradient Vv. The variation of energy dissipation with scale, which is related to 
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intermittency, will be, therefore, kept finite throughout. Then, the explicit form of jturb ma Y 
be deduced from the following dynamical requirements. 

(i) Temporal homogeneity: For turbulence driven by statistically steady forcing, the 
dynamics must be invariant under translations in time. This rules out an explicit time 
dependence of jturb- 

(ii) Spatial homogeneity: The dynamics must be invariant under translations in space. 
This excludes explicit x-dependence of V.jt ur (,. 

(iii) Isotropy: The spatial derivatives of Vt> must appear isotropically. 
We can, thus, write the following general form: 

jturb = Ai(-V)(W) 2 + \ 2 (-V) 2 (Vv) + ... , (3) 

where A's are positive constants. (The space (time) independence of A's follows from the 
invariance of V.jturb under space (time) translations, while their positivity is attributed 
to the smoothing effect of ] t urb on the spatial variations of Vi>.) Note that, because of 
homogeneity, each term in (J3J) vanishes when averaged in compliance with < j t urb >= 0, as 
required by (J2J). 

We, thus, propose the following equation as our model of steadily forced isotropic turbu- 
lence: 

d t v = vV 2 v + \ 1 W 2 {Wv) 2 -\ 2 V i v + r ] . (4) 

This equation contains all the relevant terms describing the phenomenon. Scaling arguments 
indicate that terms arising from other first and second order derivatives of Vi> in j t urb (such 
as (— V)(Vf) 4 , (— V) 2 (Vf) 3 , ...), as well as those arising from higher order derivatives, are 
all irrelevant to the scaling behavior of the equation and are, hence, disregardable. 

Let us now introduce the size, I, of turbulent eddies, which is, by definition, the length 
scale over which the velocity field varies appreciably. The lower bound to I is determined 
by the viscous dissipation length, l d , and the upper bound by the turbulent system size, L, 
which is set by the stirring forces. The competition between the linear terms V 2 t> and V 4 t> 
in P|, which represents the competition between viscous and inertial processes, generates a 



characteristic length scale l c = v /(A 2 /i/) [2l|. For I <C l c , i.e. z/|V 2 t>| A 2 |V 4 f|, turbulent 
fluctuations dominate and the dynamics of the velocity profile can be represented by ne- 
glecting viscosity. This 'small-scale' regime, thus, corresponds to the inertial regime of the 
steadily forced isotropic turbulence. The dynamics in this regime is, whence, governed by 



the nonlinear equation 

d tV = \ 1 V 2 (Vv) 2 -\ 2 V a v + t]. (5) 
(In the 'large-scale' regime, I 3> l c , viscosity dominates and the dynamics will be, therefore, 



described by the linear Edwards- Wilkinson 



22j equation. This regime is not relevant to 



inertial range dynamics.) If l c is large enough- comparable with the system size L- the 
inertial regime will involve the whole range Id < / <C L. This is generally expected to 
occur for large-scale/Reynolds-number flows. On the other hand, if l c is small, one should 
concentrate on very small eddy sizes (<C l c ) in order to observe isotropic turbulence. This, 
of course, corresponds to the local isotropy hypothesis for laboratory-scale flows, which is 
corroborated by our model. The dynamics described by equation (j5J) is independent of 
viscosity and is strongly influenced by the nonlinear term (as shown by scaling arguments). 
Viscosity plays role only in so far as to determine the extent of the inertial regime through 

lc- 



III. THE CONSEQUENCES 

The scaling behavior of the two-point correlation function < w(x, t)v (x', t) > in the 
inertial regime can be derived from a dynamical renormalization group analysis of equation 
(J5J). This equation has also made appearance in surface growth problems [l5|, where it was 
introduced by Lai and Das Sarma [2j| to describe growth by the molecular-beam-epitaxy 
(MBE) process. As mentioned in the Introduction, the turbulent velocity profile may be 
viewed as an evolving (three-dimensional) hypersurface that attains a self-affine form due 
to the roughness generated by the collective action of turbulent fluctuations. The two-point 
velocity correlation function yields a measure of roughness of this self-affine hypersurface. 
The higher order correlation functions of the model were argued by Krug to display 
multiscaling (intermittency), the origin of which was traced back to the nonlinear term. This 
reconciles with the general understanding that intermittency sets in for the fluid because of 
the nonlinearity of the Navier-Stokes equation (e.g. Q, O ur m odel of steadily forced 

isotropic turbulence, thus, falls into the same universality class of scaling behavior as the 
MBE growth (albeit) on a (non-physical) three-dimensional substrate. We have the scaling 
form 

< v(0, t)v{x, t) >= x 2a f(tx~ z ), (6) 
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where x = |x|, and the scaling function f(y) satisfies f(y ^> 1) ~ const., f(y < 1) ~ 
y 2a l z . In the language of interface problems, a and z are the roughness and the dynamics 
exponents, respectively, which are calculated [23] to be a = 1/3, z = 11/3 for dimension 
three. According to (jUJ), the velocity correlation begins to increase with time as ~ t 2a ^ z ( = 
t 2 / 11 ) but eventually reaches the stationary stage where it scales as ~ x 2a . With a = 1/3, 
this is just the Kolmogorov's 2/3 law for the fully developed (stationary) stage of isotropic 
turbulence. Of course, because Kolmogorov's theory deals only with the stationary stage 
of the fully developed turbulence, it does not at all comment on z. This exponent, which 
characterizes the time-dependent behavior of turbulence, is a dynamic feature of our model. 
As a corollary of the scaling form (jHJ), the transient time evolution and the stationary 
st at e st a ges can be characterized, respectively, by t <C r and t t, where r ~ L z . The 
crossover time, r, thus increases indefinitely with the turbulent system size implying that 
fully developed turbulence constitutes a finite-size effect. 

To understand the physical origin of the transient behavior, let us restate the above 
results in terms of the energy spectrum instead of the velocity correlation function. We 
define the Fourier transform pairs 

vV 7 k 

d{k,t) = A= J v(x,t) e" ik - x d 3 x, 



'V 

where V is the volume occupied by the fluid. One obtains, on using the homogeneity 
property, 

< |£(k,t)| 2 >= / < u(0,t)v(x,t) > e~ ik - x d 3 x 
< v(0,t)v(x,t) >= — p J < \d{k,t)\ 2 > e ik - x d 3 k, 

having let V — » oo (as required for rigorous isotropy) after taking the averages. We thus 
find, using (JHj), that 

< \v{k,t)\ 2 >= k-( 2a+ ^g{tk z ), 

where k = |k|, and the scaling function g(y) behaves asymptotically just like f(y). The 
turbulent kinetic enrgy per unit mass of the fluid, namely, 

£ = l < " 2> =l(2^/ <|8(k - i)|2># * : - 



whence yields the energy spectrum 

£(k,t) = ^ k-( 2a+1 1g(tk*), (7) 

for the inertial range k c <C k < kd, where k c = 2n/l c and kd = 2n/ld, of course. As implied 
by (JI|), the energy per unit wave-number of mode k begins to increase with time as \t 2a ^ z ', 
eventually saturating with the stationary scaling form k~( 2a+l \ Again, with a = 1/3, this 
is just the Kolmogorov's —5/3 spectrum. The transient rise in the energy spectrum can be 
understood as follows. The steady energy input from the stirring forces continually goes 
directly into the smallest k modes that are associated with the largest-scale eddies being 
produced. Since initially only these small modes are present, energy dissipation, which is 
associated with the largest modes, does not take place. The absorbed energy, therefore, 
rises continually until it has had time to cascade to the largest modes in the fully developed 
stage. The energy build-up then stops because of dissipation, resulting in the stationary 
Kolmogorov spectrum. The transient behavior is, thus, attributed to the time lag between 
the energy injection at large scales and its dissipation at fine scales j^. The larger the 
system size, the larger the range involved in the full cascade process and, hence, the larger 
the time lag is expected to be. This (together with the absence of dimensional parameters) 
explains the scaling behavior r ~ L z for the crossover time to the fully developed turbulence. 
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